Ergoregion instability of ultra-compact astrophysical objects 
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Most of the properties of black holes can be mimicked by horizonless compact objects such as 
gravastars and boson stars. We show that these ultra-compact objects develop a strong ergoregion 
instability when rapidly spinning. Instability timescales can be of the order of 0.1 seconds to 1 
week for objects with mass M = 1 — 1O 6 M0 and angular momentum J > QAM 2 . This provides a 
strong indication that ultra-compact objects with large rotation are black holes. Explosive events 
due to ergoregion instability have a well-defined gravitational-wave signature. These events could 
be detected by next-generation gravitational-wave detectors such as Advanced LIGO or LISA. 

PACS numbers: 04.40.Dg,04.30.Nk,04.25.Nx,95.85.Sz,04.70.-s 



I. INTRODUCTION 

Black holes (BHs) in Einstein-Maxwell theory are char- 
acterized by three parameters [lj: Mass M, electric 
charge Q and angular momentum J = aM ^ M 2 . 
BHs are thought to be abundant objects in the Universe 
0. Their mass is estimated to vary between 3M Q and 
10 9 5 M© or higher. They are likely to be electrically neu- 
tral because of the effect of surrounding plasma Q and 
their angular momentum is expected to be close to the ex- 
tremal limit because of accretion and merger events [H Q . 
An example of astrophysical BH is the compact primary 
of the binary X-ray source GRS 1915 + 105, which recent 
observations identify as a rapidly-rotating object of spin 
a > 0.98 M Many of the supermassive BHs which 
are thought to power quasars seem to be rotating near 
the Kerr bound [7[. 

Despite the wealth of circumstantial evidence, there is 
no definite observational proof of the existence of astro- 
physical BHs. (A review and a critique of current evi- 
dence can be found in Ref. Q and Ref. [H, respectively. 
See also Ref. [9j for a stimulating mini- review.) Astro- 
physical objects without event horizon, yet observation- 
ally indistinguishable from BHs, cannot be excluded a 
priori. 
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Dark energy stars or "gravastars" are compact objects 
with de Sitter interior and Schwarzschild exterior [Tol. [TTIl . 
These two regions are glued together around the would- 
be horizon by an ultra-stiff thin shell. In this model, a 
gravitationally collapsing star undergoes a phase tran- 
sition that prevents further collapse. The thickness of 
the shell sets an upper limit to the mass of the gravas- 
tar [ll|, [T^, [H[ . (A thorough analysis of the maximum 
compactness of gravastars can be found in Ref. [14j|.) 
Generalizations of the original model use a Born-Infcld 
phantom field [l5| , dark energy equation of state [l6| or 
non-linear electrodynamics jl7j |. Models without shells 
or discontinuities have been investigated in Ref. [HI]. 

Boson stars are macroscopic quantum states which are 
prevented from undergoing complete gravitational col- 
lapse by Heisenberg uncertainty principle [l9|, [20( . Their 
models differ in the scalar self-interaction potential plj 
and can be divided in three classes 1221: 



Miniboson stars. If the scalar field is non-interacting, the 
maximum boson star mass is M max ~ 0.633mp lanck /m 
fl9| . This value is much smaller than the Chandrasckhar 

k /m 2 . Stability of 
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mass for fermion stars, Mch 
supermassive objects requires an ultralight boson of mass 
m = 8.45 x 10" 26 GeV (10 6 M o /M max ). 

Massive boson stars. The requirement of ultralight 
bosons can be lifted if the scalar field possesses a quartic 
self-interaction potential of the form A|<^>| 4 /4 [23|. As 
long as the coupling constant A is much larger than 
(m/mpianck) 2 , the maximum boson star mass can be of 
the order of the Chandrasekhar mass or larger, M max ~ 
0.062A 1 / 2 mp lanck /m 2 . Thus supermassive objects may 
exist. Boson mass and coupling constant are related by 
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772 = 3.2 x 1(T 4 GeV A 1 / 4 (lO 6 M /M roax ) 1 /2. 

Nontopological soliton stars. If the self-interaction takes 
the form U = m 2 \4>\ 2 (l — |(/>| 2 /0q) 2 , compact nondisper- 
sive solutions with a finite mass may exist even in the 
absence of gravity [24|. The critical mass of these ob- 
jects is Af max ~ 0.0198m| lanck / ( tti4>q). If </>o ~ 777, a star 
of mass M ~ 10 6 M Q corresponds to a heavy boson of 
mass to ~ 500 GeV. 

Boson stars are indistinguishable from BHs in the New- 
tonian regime. Since they are very compact, deviations 
in the properties of orbiting objects occur close to the 
Schwarzschild radius and are not easily detectable elec- 
tromagnetically [25J, |26| . If the scalar field interacts only 
gravitationally with matter, compact objects may safely 
inspiral "inside" the boson star, the only difference with 
a BH being the absence of an event horizon [13] • Lack of 
strong constraints on boson masses makes these models 
difficult to rule out. 

Gravastars and boson stars provide viable alternatives 
to astrophysical BHs. To ascertain the true nature of 
ultra-compact objects it is thus important to devise ob- 
servational tests to distinguish these ultra-compact ob- 
jects from ordinary BHs. The traditional way to distin- 
guish a BH from a neutron star is to measure its mass. If 
the latter is larger than the Chandrasekhar limit, the ob- 
ject is believed to be a BH. However, this method cannot 
be used for the ultra-compact objects discussed above, 
because of their broad mass spectrum. A possibility is 
to look for observables related to the accretion mecha- 
nism. For example, the luminosity of quiescent BHs is 
lower than the maximum luminosity which is allowed by 
the gas present in their environment 28]. If the BH ac- 
cretion rate is much smaller than the Eddington rate, 
the radiative efficiency is also very small 29]. Another 
possibility is to exploit the absence of a boundary layer 
at the surface. Compact stars with accretion disks have 
typically a narrow viscous boundary layer near their sur- 
face, which allows the release of a considerable amount 
of heat energy. On the other hand, if the central object 
is a BH, no boundary layer is formed. Arguments of this 
kind have already excluded many gravastar candidates 
[30| . Absence of type I X-ray bursts is another powerful 
indicator of the presence of a BH. Several studies on type 
I bursts show that they are produced when gas accretes 
on the surface of a NS [31| , which then undergoes a semi- 
regular series of thermonuclear explosions. Since BHs do 
not have surfaces, the surrounding gas cannot accumu- 
late and thermonuclear instabilities do not develop. 

Another very promising observational method to probe 
the structure of ultra-compact objects is gravitational 
wave astronomy (32J. Gravitational wave detectors such 
as LIGO H|, VIRGO Q, TAMA [35], or LISA [H 
could provide an efficient way to study these objects with- 
out intervening effects due to the interstellar medium. 
For example, the inspiral process of two compact objects 
allows a precise determination of their mass [371 ] and mul- 
tipole moments [13, [H, [H, H(| ■ The gravitational wave- 



form in the presence of a surface is also expected to be 
different than the waveform in the presence of an event 
horizon [4l[. A first study on the distinctive features of 
the inspiral signal of boson stars can be found in Ref. 
[42| . Detection of gravitational resonant modes due to 
the gravitational potential well could also provide a test 
for the presence of a horizon Preliminary studies 

for gravastars indicate that this method may be very ef- 
ficient if the source is not too far away and gravitational 
wave production is significant [H, [H, [44| . 

In this paper, we propose a new method for discrim- 
inating BHs from ultra-compact horizonless objects and 
apply it to gravastars and bosons stars. Our method 
uses the fact that compact rotating objects without event 
horizon are unstable when an ergoregion is present. The 
origin of this ergoregion instability can be traced back to 
superradiant scattering. In a scattering process, super- 
radiance occurs when scattered waves have amplitudes 
larger than incident waves. This leads to extraction of 
energy from the scattering body [45], [4(| H3]. Instabil- 
ity may arise whenever this process is allowed to repeat 
itself ad infinitum. This happens, for example, when a 
BH is surrounded by a "mirror" that scatters the super- 
radiant wave back to the horizon, amplifying it at each 
scattering. The total extracted energy grows exponen- 
tially with time until the radiation pressure destroys the 
mirror in a process called BH bomb (see Refs. (48|, |49|). 
If the mirror is inside the ergoregion, superradiance may 
lead to an inverted BH bomb. Some superradiant waves 
escape to infinity carrying positive energy, causing the 
energy inside the ergoregion to decrease and eventually 
generating an instability. This may occur for any rotat- 
ing star with an ergoregion: The mirror can be either its 
surface or, for a star made of matter non-interacting with 
the wave, its center. BHs are stable, which could be due 
to the absorption by the event horizon being larger than 
superradiant amplification. 

The ergoregion instability appears in any system with 
ergoregions and no horizons [50| . (See also Ref. [HTj] for 
a exhaustive discussion.) Explicit computations for ordi- 
nary rotating stars can be found in Ref. [52l . [53[ , where 
typical instability timescales are shown to be larger than 
the Hubble time. In this case, the ergoregion instabil- 
ity is too weak to produce any effect on the evolution of 
the star. This conclusion changes drastically for ultra- 
compact stars. For compactness M > 0.5R and an- 
gular momentum J > 0.4A/ 2 , we find that instability 
timescales range approximately from 0.1 seconds to 1 
week for objects with mass in the range M ~ IMq to 
10 6 M Q , further decreasing for larger rotation rates. 

Due to the difficulty of handling gravitational pertur- 
bations for rotating objects, the calculations below are 
mostly restricted to scalar perturbations. However, we 
are able to show that the equation for axial gravitational 
perturbations of gravastars is identical to the equation for 
scalar perturbations in the large I = m limit. There are 
also generic arguments suggesting that the timescale of 
gravitational perturbations is smaller than the timescale 
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of scalar perturbations for low m. Thus our investigation 
seems to rule out some of these ultra-compact, rapidly 
spinning objects as BH candidates. 

This paper is organized as follows. In Section [ill we 
review the main characteristics of the ultra-compact ob- 
jects discussed above. Our discussion is non-exhaustive 
and strictly limited to concepts and tools which will be 
needed in the rest of the paper. Section ["II Al introduces 
the two gravastar models which will be discussed in the 
subsequent analysis. Since there are no known solutions 
describing rotating gravastars, the formalism of Refs. 

[H[ will be used to discuss rotating gravastars. Sec- 
tion III Bl introduces boson stars [H| . Numerical results 
for rotating boson stars are taken from Ref. [561 ]. Sec- 
tion Hill presents a detailed investigation of the instability 
of boson stars and gravastars using the WKB approxi- 
mation. The WKB analysis is then compared with full 
numerical results obtained by direct integration of the 
Klein-Gordon equation. Detectability of the ergoregion 
instability by gravitational-wave detectors is addressed 
in Sect. HVl Section [V] contains a brief discussion of the 
results and concludes the paper. 

Geometrized units (G = c = 1) are used through- 
out the paper, except when numerical results for ro- 
tating boson stars from Ref. [56| are discussed (Section 
IBB]) . In this case, the Newton constant is defined as 
G = 0.05/(4tt). 



consists of three regions: 

I. Interior: < r < r\ , p = —p , 

II. Shell: n < r < r 2 , p = p , (2.2) 

III. Exterior: r 2 < r , p = p = , 

where p is the energy density and p is the isotropic pres- 
sure of the gravastar. In region I, p = 3Hq/8tt is constant 
and the metric is dc Sitter: 



< r < T-i , 



(2.3) 



where C is an integration constant to be determined from 
matching conditions. In region III the spacetime is de- 
scribed by the Schwarzschild metric, 
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(2.4) 



In region II, the metric is determined by the system of 
equations, 



din r 
din h 



dh 



1 — w — h 

1 — w — h 
l + w-3h 



1 

B 



d\n w . 



(2.5) 
(2.6) 



II. STRUCTURE OF ULTRA-COMPACT 
ASTROPHYSICAL OBJECTS 



This section discusses the main properties of gravastars 
and boson stars. The derivation of nonrotating solutions 
is partly based on Refs. [ll|, [l3|, [23|, [56j . 



A. Gravastars 

Although exact solutions for spinning gravastars are 
not known, they can be studied in the limit of slow rota- 
tion by perturbing the nonrotating solutions [Hi] ]. This 
procedure was used in Ref. [55j to study the existence of 
ergoregions for ordinary rotating stars with uniform den- 
sity. Their analysis is repeated below for gravastars. In 
the following, we discuss the original thin-shell model by 
Mazur and Mottola [ll| and the anisotropic fluid model 
by Chirenti and Rezzolla [H, [l8| . 



where w = 8irr 2 p and wf /r 2 = constant. A simple ana- 
lytical solution can be obtained for a thin shell [llj. In 
the limit n — > r 2 , one obtains 



1 (1 + u;) 2 , 

- ~ 6± '- « 1 

B w 



(2.7) 



where e is an integration constant. The continuity of 
the metric coefficients / and B at r\ and r 2 implies that 
e, C, M and Hq are related to r%, r 2 , w\ = w(r\) and 
w 2 = w(r 2 ) by [131 



_i n ra(i n !^_J_ 

Wl w 2 



C 
M 



n 

1 + w 2 

1 + W\ 



1 

U>1 



[2 

2 



1 - 



1 - 



e{l + w 2 f 
w 2 

e{l + w±f 



Wl 



(2.8) 
(2.9) 
(2-10) 
(2.11) 



1. Nonrotating thin-shell model 

In this model, the spacetime 

ds 2 = -f(r)dt 2 + B{r)dr 2 + r 2 dVi\ (2.1) 



The above relations and Eq. (|2.5p completely determine 
the structure of the gravastar. A typical solution is shown 
in the upper panel of Fig. [T]for r 2 = 1.05, r*i = 1, Wi = 
350 and w 2 = 1. 
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FIG. 1: Top panel: Metric coefficients for the thin-shell model 
gravastar with V2 = 1.05, n = 1, wi = 350 and W2 = 1, corre- 
sponding to M ~ 0.485 ri. Bottom panel: Anisotropic pres- 
sure model for V2 = 2.2, ri = 1.8 and M = 1, corresponding 
to A'/ - 0.45r 2 . 



2. Nonrotating gravastars with anisotropic pressure 

This model assumes a thick shell with continuous pro- 
file of anisotropic pressure to avoid the introduction of 
an infinitesimally thin shell. The stress-energy tensor is 
T^ v — diag[— p,p r ,Pt,Pt]i where p r and p t are the ra- 
dial and tangential pressures, respectively. The density 
function is 

po , < r < ri region I 

p(r) = I ar 3 + br 2 + cr + d , r\ < r < r 2 region II 
, r 2 < r region III 

with boundary conditions p(0) = p{r\) = po, p{i~2) = 
p'(ri) = p'{r 2 ) = and 



2po 



(r 2 - ri) 3 ' 
6porir 2 
(r 2 - ri) 3 ' 



/| = _ 3p (r 2 + n) t (2i2) 



(r 2 -ri) 3 
p (r 3 - 3rir|) 
(r 2 - ri) 3 



(2.13) 



The density is related to the total mass M by 

Po_ 15 

M ~ 2Tr(ri + r 2 ){2r\ + r x r 2 + 2rf) ' 

The radial pressure p r is chosen as [l3| 

2^ 



a — (1 + a) 



(2.14) 



(2.15) 



where the parameter a is determined by demanding that 
the maximum sound speed coincides with the speed of 
light. (This requirement rules out superluminal behavior 
and implies a ~ 2.21.) The metric coefficients are 



f=(l- — ) e r «- r (-> , I = 1 - 



(2.16) 



where 



and 



m(r) = / A-nr 2 pdr 
Jo 



r(r) 



2m(r) + 87rr 3 p T 
r(r — 2m(r)) 



(2.17) 



The above equations completely determine the structure 
of the gravastar. Both the metric and its derivatives are 
continuous across r 2 and throughout the spacetime. The 
behaviors of the metric coefficients for a typical gravastar 
are shown in the bottom panel of Fig. [TJ 



3. Slowly rotating rigid gravastars and ergoregions 

There are no known solutions describing rotating 
gravastars. Thus an analysis of the ergoregion instabil- 
ity for these objects is nontrivial. Fortunately, slowly 
rotating solutions can be obtained using the formalism 
developed in Ref. [54|, which we now extend to the case 
of anisotropic stresses. 

A rotation of order fl gives corrections of order Q 2 in 
the diagonal coefficients of the metric (|2.1[) and intro- 
duces a non-diagonal term of order O, 



9t</> = -Qg<t>4> ■ 



(2.18) 



where <p is the azimuthal coordinate. The metric coeffi- 
cient gtcf, defines the angular velocity of frame dragging 
C = C( r )- The full metric is 



ds 2 = -f{r)dt 2 + B{r)dr 



l d0 2 



+ r 2 sin 2 9 (# - C{r)dt) 2 . (2.19) 
We consider the anisotropic fluid stress-energy tensor 

= (p + PtWU v + ptgT + (Pr ~ Pt)s»s v , (2.20) 
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where 



jjr 

u* 



„-0, 



U e = 0, U' t ' = W t , 

[- {gtt + 205^ + fi 2 ^)] 



-1/2 



Equation (|2.20p describes an anisotropic fluid with radial 
pressure p r and tangential pressure rotating with an- 
gular velocity £1 as measured by an observer at rest in the 
(t , r , 6 , 4>) coordinates. If the gravastar rotates rigidly, 
i.e. = constant, the Einstein equations at order £ give 



- 8np 
8irp r 

8irp t 



B - B 



rB' 



f-Bf + rf> 



r 2 Bf ' 
2f 2 B' + rBf 2 

ArB 2 f 2 
f(rB'f'-2B(f' + rf")) 
ArB 2 f 2 

An equation for Q(f) is obtained by considering 
R t( p — 8tt \ Tt<j> — -^9t<$>T 



(2.21) 
(2.22) 

(2.23) 

(2.24) 



Using Eqs. (j2~2T|) and (f2~22|) . Eq. (|2~24|) is written as 

C" + C'(J: + j) = 167rB(r)(C-n)(p + pt) , (2.25) 

where j = (fB)^ 1 / 2 is evaluated at zeroth order and 
p, pt are given in terms of the nonrotating geometry by 
Eq. (|2.21|) and Eq. (|2.23p . respectively. The above equa- 
tion reduces to the corresponding equation in Ref. [54| for 
isotropic fluids. Solutions of Eq. (|2.25[) describe rotating 
gravastars to first order in Q. 

Spinning gravastars may possess ergoregions. A sim- 
ple but general procedure to determine their presence 
for slowly rotating stars is described in Ref. [55]. This 
method requires only a knowledge of the metric of non- 
rotating objects and compares favorably with more so- 
phisticated numerical analyses [57l ]. 

The ergoregion can be found by computing the surface 
on which g tt vanishes (55| : 



= -f(r) + ( 2 r 2 sin 2 



(2.26) 



Equation (|2.26[) is expected to be a good approximation 
to the location of the ergoregion specially for very com- 
pact stars UK . The solution of Eq. (|2.26[) is topologically 
a torus. In the equatorial plane we have 



(2.27) 



The existence and the boundaries of the ergoregions can 
be computed from the above equations. Equation (|2.25p 
is integrated from the origin with initial conditions (fl — 



CY = and (f2 — £) finite. Changing the value of (fi — 
(), the whole space of slowly rotating gravastars can be 
obtained. The exterior solution satisfies — £ = Sl(l — 
21 /r 3 ), where / is the moment of inertia of the star. 
Demanding the continuity of both (17 — £)' and (fi — £), C 
and / are uniquely determined. The rotation parameter 
n depends on the initial condition at the origin. 





r/r, 



FIG. 2: Top panel: Size of the ergoregion for three differ- 
ent gravastars in the thin-shell limit. The vertical axis gives 
the angular momentum of the gravastar in units of its total 
mass. The horizontal axis gives the locations of the ergore- 
gion boundaries in units of the gravastar radius T2- Each curve 
refers to a different gravastar. The minima of the curves de- 
termine the existence and extent of ergoregions. The size of 
an ergoregion can be found by drawing a horizontal line at a 
given value of J/M 2 ; its intersections with the curves give the 
radii of the ergoregion boundaries. From top to bottom the 
three curves refer to V2 = 1.3, r\ = 1, ioi = 50 and W2 = 1, 
corresponding to M ~ 0.39 V2\ T2 = 1.2, ri = 1, w\ = 150 
and W2 = 1, corresponding to M ~ 0.44 T2\ r2 = 1.05, r\ = 1, 
wi — 350 and W2 = 1, corresponding to M ~ 0.49 T2- Bot- 
tom panel: J/M 2 and angular frequency Q, for the anisotropic 
pressure model with T2 = 2.2, r\ = 1.8 and M — 1. The an- 
gular frequency is always very small, thus the slow-rotation 
formalism applies. These results extend up to the Keplerian 
frequency Qk- 



6 



Figure [2] shows the results for three different gravas- 
tars described in the previous sections. The minima of 
the curves are the minimum values of J/M 2 which are 
required for the existence of the ergoregion. Comparison 
with the results for stars of uniform density [55j , shows 
that ergoregions form more easily around gravastars due 
to their higher compactness. Figure [2] also shows that 
the ergoregions spread inside the gravastar. (The ergore- 
gion can be located by drawing an horizontal line at the 
desired value of J/M 2 , as explained in the caption.) 

Gravastars spinning above a given threshold are not 
stable against mass shedding 58] . Instability arises when 
the centrifugal force is strong enough to disrupt the 
star. In Newtonian gravity, the equatorial mass shed- 
ding frequency is approximately the Keplerian frequency 
MVL K = (M/i?) 3 / 2 . Although corrections to the Ke- 
plerian frequency are expected in a general relativistic 
framework, provides a good estimator for the validity 
of the slow-rotation approximation. (See Ref. [59| for a 
comparison of the slow-rotation regime vs. full numerical 
results.) In the following, the slow-rotation approxima- 
tion will be considered valid for VI/VLk < 1- Numerical 
results extend up to Q ~ Qk- 



B. Boson stars 

A well-known example of nonrotating boson star is the 
model by Colpi, Shapiro and Wasserman (CSW) [23j ] - 
A variation of the CSW model which allows for rotat- 
ing solutions is the Kleihaus, Kunz, List and Schaffer 
(KKLS) model 0. The KKLS solution is based on the 
self-interacting complex scalar field $ with Lagrangian 
density 

C K KLS = -\g^{®*^,u + -U(\$\), (2.28) 

where U(\$\) = A|$| 2 (|$| 4 - a|$| 2 + b). The mass of 
the boson is given by ma = VXb. The equations for the 
boson star structure can be solved by setting 



ds 2 = -fdt 2 + - 



g (di 



2 , „2 dQ 2 



f 

+ r 2 sin 2 9 (dip - C(r) dtf 



(2.29) 



and <& = 4> giust+incp^ where the me tric components and 
the real function <f> depend only on r and 9. The require- 
ment that $ is single- valued implies n — 0, ±1, ±2, 

The solution has spherical symmetry for n = and axial 
symmetry otherwise. The mass M and the angular mo- 
mentum J can be read off from the asymptotic expansion 
of / and (, 

M=^- lim r 2 d r f , J = -L lim r 3 C , (2.30) 

l(j i — >oo ZLr r^oa 

respectively. Since the Lagrangian density is invariant 
under a global £7(1) transformation, the current — 



-i$*<9 Al< l> + c.c. is conserved. The associated charge is 

Q = Airu s [°° f | 5 | 1/2 i + — -Vdr<». (2.31) 
Jo Jo J V u s r J 



It can be shown that the angular momentum J and the 
scalar charge Q of stationary solutions, i.e., solutions 
with a timelike and rotational Killing vector, are related 
by the quantization condition J = nQ [60| . This may 
erroneously lead to the conclusion that boson stars are 
gravitationally stable because a continuous extraction of 
angular momentum is needed to trigger instability. How- 
ever, the above quantization condition applies for objects 
for which such Killing vectors can be defined and we will 
assume that this condition is broken in the presence of 
perturbations. In this case, an arbitrary amount of an- 
gular momentum can be extracted and it will be shown 
below that gravitational- wave emission leads to an insta- 
bility on very short timescales. 

The numerical procedure to extract the metric and 
the scalar field is described in Ref. [H(|. Through- 
out the paper we will consider solutions with n = 2, 
b = 1.1, A = 1.0, a — 2.0 and different values of 
(J,M) = (3781, 1296), (3400, 1081), (2800,906), corre- 
sponding to J/{GM 2 ) ~ 0.566, 0.731 and 0.858, respec- 
tively. The n = 1 solutions in Ref. [56l | exhibit similar 
features. The two top panels of Fig. [3] show the met- 
ric functions for boson stars with J/(GM 2 ) ~ 0.566 and 
0.858 along the equatorial plane. The change in the met- 
ric potentials from 9 = tt/2 to 9 = ir/4 for these so- 
lutions is plotted in the bottom panels of Fig. [3] The 
metric functions do not depend significantly on the lon- 
gitudinal angle. Figure Q] gives g u as a function of dis- 
tance for the case with J/(GM 2 ) ~ 0.566 at the equa- 
tor. The behavior of gu demonstates that boson stars 
develop ergoregions deeply inside the star. For this par- 
ticular choice of parameters, the ergoregion extends from 
r/(GM) ~ 0.0471 to 0.770. A more complete discussion 
on the ergoregions of rotating boson stars can be found 
in Ref. M- 



III. ERGOREGION INSTABILITY FOR 
ROTATING STARS 

The stability of ultra-compact objects can be stud- 
ied perturbatively by considering small deviations around 
equilibrium. As explained in the introduction, we con- 
sider only scalar perturbations. This is justified as fol- 
lows. Axial gravitational perturbations are described in 
the large I = m limit by the same equation of scalar 
perturbations. In this regime, our results describe both 
kinds of perturbations. In the small I = m limit, gravita- 
tional perturbations are expected to have shorter growth 
times than scalar perturbations: black holes are charac- 
terized by a superradiant amplification of spin-2 fields 
which is much stronger than the superradiant amplifica- 
tion of other fields. This is due to the potential barrier 
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FIG. 3: Left panel: Metric coefficients for a rotating boson star along the equatorial plane, with parameters n = 2, b = f .f , 
A = f.O, a = 2.0, J/{GM 2 ) ~ 0.566. Right panel: Fractional difference of the metric potentials between 9 = 7r/2 and 6 = 7r/4 
for the same star. The plot gives the maximum possible fractional difference between these quantities. 
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FIG. 4: The gtt metric coefficient for a boson star with 



J/{GM 2 



0.566 at its equator. The ergoregion is identi- 



fied by the region inside the dotted vertical lines and extends 
from r/(GM) ~ 0.047 to 0.770. 



outside the horizon having different behavior for differ- 
ent spin-field s. Since the ultra-compact objects we are 
dealing with are also characterized by a relativistic poten- 
tial barrier, gravitational perturbations are expected to 
couple more strongly to the ergoregion and have smaller 
instability timescales. This conclusion is also verified un- 
der certain simplifying assumptions in Ref. [6l| . Thus, 
scalar perturbations should provide a lower bound on the 
strength of the instability. 



A. Axial gravitational perturbations for perfect 
fluid stars 

In the large I = m regime, axial gravitational per- 
turbations [6l|, are described by a simple equation. 
In general, axial and polar perturbations are coupled 
when rotation is included [63j |. For simplicity, we will 
assume that the zeroth-order polar perturbations vanish 
and there is no coupling. The full metric is a perturba- 
tion of Eq. (|2T9|) 



ds 2 = dsn 



2j2(h l o m (t,r) + h[ m (t,r)) 



(- 



,-1 



YimdS + smOdaY, 



,(3.1) 



where ds^ is the unperturbed metric (|2.19[) and Yj m are 
scalar spherical harmonics. The quantities 



Am 



Vf(r)/B(r)K 6 , h l { n = y/WHW) ^ 



(3.2) 

satisfy the system of equations (see Eqs. (13)-(16) in 
Ref [63) 



K' 3 = 16n(p + p)u 3 — 



2K* 



l 2 + l-2 T ^ 2m('V 4 
o K & - 



(3.3) 



B ( f 1 B f 2 N 
KL = (-uj + mC)Vi - )K 6 



BK 3 , 
(3.4) 

where K3 and Kq are two extrinsic curvature variables 
and 



"3 



I 2 + 1-2 



(-w + mC)K 3 
2m(Q - c) 



2m(Q - C) - 1(1 + + mfi) 



2m(' K Q 
l(l + l)~B 

K 6 . 



, (3.5) 
(3.6) 
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Here and throughout this paper, it is understood that 
u) = ivi m . In the large I = m limit, Eqs. (|3. 3|) and (|3.4[) 
reduce to 



2 171 

-K 3 + —K 6 , 



(3.7) 



if' = (E + OV.K3 - # 6 , 

(3.8) 

where E = —ujjm. Combining Eqs. (|3.7[) - (|3.8[) and ne- 
glecting terms of order 1/m 2 , it follows 



K'! + m 2 j ((S + C) 2 -i) ^ = 0. C-!. 1 )! 



This equation also describes the scalar perturbations of 
gravastars, as it will be shown below. 




r/r, 



B. Scalar field instability for slowly rotating 
gravastars: WKB approach 

Consider now a minimally coupled scalar field in the 
background of a gravastar. The metric of gravastars is 
given by Eq. (|2.19|) . In the large I = m limit, which is 
appropriate for a WKB analysis [H, , the scalar field 
can be expanded as 



( - im 



$ = E^)e XP [-5/(' + ^ + g) dr 

Ira 7 

(3-10) 

The functions xim are determined by the Klein-Gordon 
equation, which yields 




Xim + m 2 T(r , E)xim = , 
where E is defined as below Eq. (|3.8[) and 



T 



V± 



B(r) 
f(r) 

-c± 



(s-y+)(E-y_) , 
V7M 



(3.11) 



(3.12) 



(3.13) 



Equation (|3.11[) follows from Eq. (|3.10[) when terms of 
order 0(l/m 2 ) are dropped. Equation (|3 . can be 
shown to be identical to Eq. (|3.9[) for the axial gravita- 
tional perturbations of perfect fluid stars. Therefore, the 
following results apply to both kinds of perturbations. 

The eigenfrequencies of Eq. (13.111) c an be computed in 
the WKB approach following Ref . [52| . This method is in 
excellent agreement with full numerical results [53l . [6J] . 
The quasi-bound unstable modes are determined by 



and have timescale 



nir , n = , 1 , 2 . 



(3.14) 



r = 4 exp 



2m 



T\dr 



VTdr, (3.15) 



r/r n 



FIG. 5: Top panel: Potentials V± for the thin-shell gravastar 
with r2 = 1.3, ri = 1, w% = 50 and W2 = 1. The ergoregion 
extends from r ~ 0.247r2 to 0.832r2 and corresponds to a 
gravastar with J ~ 0.333M 2 and MQ. ~ 0.105. Bottom panel: 
Potentials for the anisotropic pressure gravastar with V2 = 
2.2, n = 1.8 and M = 1. The ergoregion extends from r ~ 
0.270r2 to r ~ 1.055r2 and rotates with angular momentum 
J/M 2 = 1.00, corresponding to Q ~ 0.250. 



where r a , are solutions of V+ = E and r c is determined 
by the condition V- — E. 

The potentials V± are displayed in Fig.[5]for the gravas- 
tar models of Sect. Ill Al The top panel shows the po- 
tential for the thin-shell model with T2 = 1.3, r*i = 1, 
w\ = 50 and W2 = 1. The gravastar rotates with 
angular frequency f2 ~ 0.105 and the ergoregion lies 
in the region r ~ (0.247, 0.832)r 2 . The bottom panel 
refers to the anisotropic pressure model with r 2 = 2.2, 
n = 1.8 and M = 1. The ergoregion extends from 
r ~ 0.270 r 2 to r ~ 1.055 r 2 and rotates with angular 
frequency ft ~ 0.250. 
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FIG. 6: Details of the ergoregion instability (m = 1 and 
m — 4) for the thin-shell gravastar of Sect. Ill A II with 
r2 = 1.3, n = 1, Wi = 50 and u>2 = 1. The plot shows 
the logarithm of the dimensionless instability timescale t/M, 
the dimensionless angular velocity MQ and the oscillation fre- 
quency E vs. the angular momentum per unit mass, J/M 2 . 
The m — 2,3 modes follow a similar pattern and are in be- 
tween the m — 1, 4 results. 



The results of the WKB computation are shown in 
Fig. [5] and Tables Hl lIIII Figure [5] displays the results for 
the least compact thin-shell gravastar of Sect. Ill Al with 
to = 1, 4. Although the WKB approximation breaks 
down at low to values, these results still provide reliable 
estimates [52j . This claim will be verified in Sect. IIII B II 
with a full numerical integration of the Klein-Gordon 
equation. Table Q] compares three different gravastars for 
J/M 2 = 1 and m = 1 , 2 , ... 5. The results show that the 
instability timescale decreases as the star becomes more 
compact. Table [TT1 refers to the most compact thin-shell 
gravastar for various angular frequencies. The instability 
timescale depends strongly on the rotation. A fit for the 
instability timescale in powers of J/M 2 yields 



logr/Af ~ a + b^J/M 2 + cJ/M 2 



(3.16) 



where a = 68.0, b = —76.7, c = 26.2 for m — 5 and 
a = 55.8, b = —61.4, c = 21.0 for m = 4, respec- 
tively. The results for lower to values show similar be- 
haviors. For instance, for to = 1 the coefficients are 
a = 19.7,6 = -16.2 and c = 5.6. Table OH shows the 
WKB results for the anisotropic pressure model for dif- 
ferent values of J/M 2 . Larger values of J/M 2 make the 
star more unstable. The instability timescales are fitted 
by 



log t/M 



b (J/M 2 



where a = —21.9, b = 39.2, c = —0.39 for to = 5 and 
a = -13.7, b = 28.5, c = -0.43 for m = 4. The trend 

TABLE I: WKB results for the instability of rotating thin- 
shell gravastars with J/M 2 — 1. The instability grows with 
compactness. 



(-10 2 ME, t/M) 



r 2 = 1.3 
9.2,6.85 x 10 3 
11,1.23 x 10 5 
2.31 x 10 6 
4.34 x 10 7 



12. 
12. 
12. 



r 2 = 1.2 
17,5.26 x 10 3 
18 , 5.63 x 10 4 
18, ,6.06 x 10 5 
18,6.58 x 10 6 



8.26 x 10 s 18,7.13 x 10' 



r 2 = 1.05 
23,9.16 x 10 3 
23,7.37 x 10 4 

23.5.98 x 10 5 
23 , 4.86 x 10 6 

23.3.99 x 10 7 



for lower-TO is similar. The coefficients for to 
a = 7.4, b = 0.57, c= -3.0. 



1 are 



Both models have similar low-m behaviors. It will be 
shown in Sect. IIII B II that the WKB results for the insta- 
bility timescale differ from the numerical results by about 
one order of magnitude at low to. On the contrary, the 
resonant frequencies match well the WKB results even 
for 1ow-to modes. Calculations show that the resonant 
frequency is Re(w) ~ aQ, where a ~ 1.1 - 1.2. 

The maximum growth time of the instability is of the 
order of a few thousand M, at least for large J. This 
instability is crucial for the star evolution. Comparison 
of Table|V]with Tables 1 and 2 of Ref. 52] shows that the 
ergoregion instability of gravastars is stronger than the 
ergoregion instability of uniform density stars by many 
orders of magnitude. This seems to be a general feature 
of all ultra-compact objects discussed here. Gravitational 
perturbations are expected to be even more unstable. 



Comparison with numerical results 



Accurate computations of the instability require nu- 
merical solutions of the Klein-Gordon equation. How- 
ever, the WKB approximation provides reliable esti- 
mates of the numerical results. The exact potential T 
of Eq. (j3~T5j) is 



T 



l(l + l)B B(u-mCY 



5B 



16B 2 



2rf 



f 

B'f 

8Bf ■ 16/ 



+ 



3/' : 



b' ir_ 

27b + 4B + 
f" 



(3.17) 



The results of the numerical integration for the 
anisotropic pressure gravastar with r-i — 2.2, r\ = 1.8, 
M = 1 and J/M 2 = 1 are shown in Table HVl 
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TABLE II: WKB results for the instability of r otatin g thin-shell gravastars with r-2 = 1.05 and ri = 1.0. The mass shedding 
limit corresponds to fi/fijr = 1, where Qk = \J M/r\. 



J/M 1 = 0.40 
Q/Q.K = 0.28 
8.4, 1.19 x 10 5 
8.7, 1.13 x 10 7 
8.8,1.08 x 10 9 
8.9, 1.04 x 10 11 



J/M' 2 = 0.50 
SI/VLk = 0.35 
10,6.28 x 10 4 
11,3.24 x 10 6 
11, 1.68 x 10 s 
11,8.75 x 10 9 



0.60 
= 0.42 



J/M' z 

13,3.70 x 10 4 
14, 1.15 x 10 6 
14,3.59 x 10 7 
14, 1.13 x 10 9 



{-\Q 2 MT,,t/M) 



J/M' 2 

16,2.38 x 10 4 
16,4.82 x 10 5 
16, 9.86 x 10 6 
16,2.03 x 10 s 



0.80 J/M' 2 
0.57 n/n K 



8.9,1.01 x 10" 11,4.58 x 10 11 14,3.59 x 10 lu 16,4.16 x 10 



0.70 TJW 

0.50 n/n* 

18, 1.64 x 10 4 20 
18,2.31 x 10 5 21 
18,3.30 x 10 6 21 
18,4.69 x 10 7 21 
18,6.75 x 10 8 21 



0.90 J/M 2 = 1.0 
0.64 n/fix = 0.71 
1.19 x 10 4 23,9.15 x 10 3 
1.24 x 10 5 23, 7.38 x 10 4 
1.30 x 10 6 23,5.97 x 10 5 
1.37 x 10 7 23,4.86 x 10 6 
1.45 x 10 s 23,3.98 x 10 7 



TABLE III: WKB results for the instability of rotating anisotropic pressure gravastars with r-z = 2.2, r\ = 1.8 and M = 1. 



(-W' 2 MT,,t/M) 



J/M 2 = 0.40 J/M 2 = 0.50 J/M 2 



Q/Q, K = 0.33 
0.24, 1.33 x 10 7 

3.1,8.25 x 10 7 
4.2, 1.31 x 10 10 
4.7,2.50 x 10 12 
5.1,5.06 x 10 14 



0.60 

Q./Q.K = 0.41 n/Q K = 0.49 
2.7, 1.13 x 10 s 5.2,2.78 x 10 4 
5.6,6.20 x 10 6 8.1, 1.14 x 10 6 



J/M r = 0.70 
n/Q K = 0.57 
7.6,1.15 x 10 4 
10,3.13 x 10 5 



6.6, 5.44 x 10 s 
7.2,5.13 x 10 
7.6,4.99 x 10 1 



in 



J/M 2 = 0.80 J/M 2 = 0.90 J/M 2 = 1.0 
n/Q K = 0.65 0,/Q.k = 0.74 {l/Q K = 0.82 



10,5.99 x 10 J 
13, 1.11 x 10 5 



9.1,5.65 x 10' 
9.7,2.95 x 10 9 
10, 1.59 x 10 11 



12,9.40 x 10° 14,2.25 x 10° 



12,3.58 x IO 3 

15.4.81 x 10 4 

17.6.82 x 10 5 



12,3.10 x 10 8 
13,9.82 x 10 9 



15,4.81 x 10 7 17, 1.02 x 10 7 
15,1.02 x 10 9 17,1.52 x 10 8 



15,2.34 x 10 d 
18,2.33 x 10 4 
19,2.45 x 10 s 
20,2.73 x 10 6 
20,3.07 x 10 7 



The WKB approximation for the real part of the fre- 
quency shows a remarkably good agreement with the 
numerical results even at low values of to. For any 
to > 1 this agreement is better than 5%. The insta- 
bility timescale seems to be more sensitive to the details 
of the WKB integration, with a level of agreement simi- 
lar to that reported in Ref. [53|]. The above results show 
that the WKB approximation correctly estimates the in- 
stability timescale for all values of to within an order of 
magnitude. 



Scalar field instability for rotating boson stars: 
WKB approach 



TABLE IV: Comparison between analytical and numerical 
results for anisotropic pressure gravastars with J/M 2 — 1, 
T2 = 2.2, ri = 1.8 and M = 1. The numerical results for the 
real part are in good agreement with the WKB results. The 
agreement is better for larger values of m. The imaginary 
parts agree within an order of magnitude. 



(-10 2 Y,M ,t/M) 



Analytical (A) 



15.0,2.34 x 
18.0,2.34 x 
19.0,2.49 x 
19.6,2.74 x 
19.8,3.07 x 



nr 
io 4 
io 5 
io (i 
io 7 



Numerical (N) 
21.9,2.47 x 10 3 
18.6,1.81 x 10 5 
18.7,2.53 x 10 6 
19.0,3.33 x 10 7 
19.3,3.76 x 10 s 



31.5% 
3.2% 
1.6% 
3.2% 
2.6% 



Consider a scalar field ^ minimally coupled to the ro- 
tating boson star geometry (not to be confused with the 
scalar field which makes up the star). Since the metric co- 
efficients depend on r and 9, the Klein-Gordon equation 
cannot be reduced, in general, to a one-dimensional prob- 
lem. Separation of variables can be achieved by requiring 
.9 = 1,/ = /(r,7r/2), k = fc(r,7r/2) and C = C(r,7r/2). 
These assumptions are justified as follows. Firstly, the 
metric function g is very close to unity throughout the 
entire coordinate region, as can be seen in Figure [3] 
Since the Klein-Gordon equation does not depend on the 
derivatives of g, it seems safe to set g = 1 in the whole 
domain. Secondly, the angular dependence of the metric 
coefficients is negligible for slow rotations. The largest 
variation for one revolution around the star is that of /, 



which is less than 100% for most cases (see Fig.[3|). More- 
over, this dependence is extremely weak for most of the 
values of r. Thirdly, perturbations are localized around 
the equator in the large I = m behavior. Therefore, eval- 
uating the metric coefficients at the equator may provide 
a reasonable approximation to the problem. The results 
obtained in this section are expected to give an order of 
magnitude estimate for the instability. 

The equation for the scalar field is obtained by expand- 
ing the latter as 



* = ^*i m (r)exp -- J 

lm 



dr 



(3.19) 
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r/(GM) 



The maximum growth time for this boson star model is 
of the order of 10 6 M for = 0.857658. This is a rela- 
tively short instability timescale, corresponding to about 
one second for a one solar mass boson star. 



IV. DETECTABILITY BY EARTH- AND 
SPACE-BASED GRAVITATIONAL WAVE 
DETECTORS 

The ergoregion instability may be of interest for grav- 
itational wave astronomy. Contrary to the r-mode insta- 
bility of neutron stars [H, [67J , the ergoregion instability 
is not limited to solar mass objects. Thus chances of de- 
tection are larger because the signal can fall in frequency 
bands where the detectors are more sensitive. 



A. Signal-to-noise ratio 



FIG. 7: Potentials V± for the boson star model with 

T 1\ 



J/(GM Z ) = 0.566. 
0.0478 to 0.779. 



The ergoregion extends from r/{GM) 



In the large I = m limit, is determined by 



0. 



where 



T = (S - V+) (S - VJ) 



V± = C± 



m 



(3.20) 

(3.21) 
(3.22) 



and terms of order O (l/m 2 ) have been neglected. The 
potentials V± are plotted in Fig. for the rotating bo- 
son star with J/(GM 2 ) = 0.566. The results of the 
WKB computation are summarized in Table [V] for m = 
1,2,... 5. An interesting feature is that the instabil- 



Detectability depends only on the energy released and 
the detector frequency bandwidth. The sky-averaged 
signal-to-noise ratio (SNR) is (68] ] 



P 



5tt 2 L> 2 



df 



1 



dE 



f 2 S h (f) df 



(4.1) 



where D is the distance to the source and Sh(f) is the 
noise power spectral density of the detector. Using dE = 
2nfdJ/m 1 the SNR for the I = m = 2 mode, i.e. the mode 
for which the instability is expected to be stronger, is 



5ttD 2 



df 



1 



dJ 



fS h (f) df 



(4.2) 



Equation (|4.2j) agrees with results in the literature [65j, 
l66l . |67| and is independent of the perturbation amplitude 
at lowest order. Higher-order corrections, however, would 
contain a dependence of the SNR ratio on the amplitude 
evolution. Fitting the resonant frequencies to 



Re[w] = 2vr/ « aO . 



(4.3) 



TABLE V: Instability for rotating boson stars with parame- 
ters n — 2, b = 1.1, A = 1.0, a = 2.0 and different values of J: 
J/(GM 2 ) = 0.566, 0.731 and 0.858. The Newton constant is 
defined as 4nG = 0.05. 



— T 

GM 



-W*GMT,,t/(GM) 



0.566139 
31,8.847 x 10 2 
36 , 7.057 x 10 3 
37,6.274 x 10 4 
38,5.824 x 10 5 

38,5.554 x 10 6 18,2.915 x 10 



= 0.730677 
6.6,6.303 x 10 3 
13,5.839 x 10 4 
16 , 9.274. x 10 5 
17,1.603 x 10 7 
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( ;m- 



0.857658 



0.68 , 1.478 x 10 b 
3.4,2.815 x 10 8 
4.9,2.815 x 10 10 
5.7,1.717 x 10 12 



ity timescale increases with the star angular momentum. 
Due to lack of sufficient numerical results we do not yet 
know whether this behavior holds for all values of the 
angular momentum or the trend changes at some point. 



one finds a ~ 1.1-1.2. From Eq. (JOJ) it follows dJ/df » 
2nl/a. Assuming the moment of inertia I ~23M 3 to be 
roughly independent of the angular velocity [6 51 ] (compu- 
tations show that /3 ~ 1 for gravastars to a very good 
approximation) Eq. (|4.2| can be rewritten as [H, \6a . [6tJ 



8/3M 3 
5aL> 2 



df 



1 



fS h (f) 



(4.4) 



The minimum and maximum frequencies in the above in- 
tegral are chosen as / m j n = 0.9/ max (we here implicitly 
assume that the amplitude of the perturbation is suf- 
ficiently large so that / m i n is reached within the typical 
operation time of a detector such as LIGO or LISA, ie ~ 
1 year). This is a conservative estimate based on a simple 
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TABLE VI: SNR for ergoregion instability of an object at 
distance of 20Mpc for LIGO/Advanced LIGO and LISA. (See 
text for details.) 



p 


M/M & 


10 


20 


30 


40 


50 


100 


LIGO 


1 


5 


11 


18 


24 


23 


Adv. LIGO 


9 


84 


193 


249 


304 


468 



P 


M/Mq 


10 4 


10 b 


10 b 


10' 


10 s 


LISA 


0.07 


21 


2428 


1469 


1417 



model for the evolution of the system. SNRs for objects 
at 20Mpc distance for LIGO/Advanced LIGO and LISA 
are shown in Table IVII for QM — 0.2. Solar-mass ob- 
jects are difficult to detect, although LIGO (Advanced 
LIGO) could be able to detect objects with M > 3OM 
(M > 10A/q). SNRs of several thousands are easily 
achieved for supermassive objects. These objects could 
be easily observed by LISA. 



B. Waveforms 

The expression for the SNR derived above is the opti- 
mal SNR. Search and detection techniques, i.e. matched 
filtering, usually require accurate theoretical templates. 
The derivation of accurate waveforms for the ergoregion 
instability is beyond the purpose of this paper. The 
physics involved is too complex and even the evolution 
of the instability itself is at present not known. The sole 
purpose of this section is to sketch the evolution of the 
most important quantities of the process. 

The instability proceeds in two steps: A phase char- 
acterized by an exponential growth, where the linear ap- 
proximation is valid, followed by a nonlinear phase. Con- 
trary to the r-mode instability, the ergoregion instabil- 
ity does not couple strongly to the fluid composing the 
object. Therefore, the nonlinear phase is expected to 
be somewhat different from the r-mode saturation phase 

[Hlzl. 

As an illustration of waveform estimation, consider the 
instability triggered by a particle in circular orbit around 
the compact object. Focusing on the I = 2 mode, the 
metric perturbations in the linear perturbation regime 
have the form 

h+ = ^h e t l T sm{ut-2<t>) X+ , (4.5) 
hx = ^h e t/T cos (tvt- 20) X x , (4.6) 
where D is the distance to the source, h << l and 
cos 2 9 + 1 



X+ 



x> 



cos 8 . 



(4.7) 



The above waveforms mimic the Newtonian waveform 
produced by a small mass in circular orbit around the 



ultra-compact object. The linear perturbation regime 
corresponds to x < 1. The exact functional form of 
the waveform (|4.7|1 is required to determine polarization 
and phase content of the waveform, but not the evolution 
of its amplitude. 

From the discussion in the previous sections, the fre- 
quency of the wave is w ~ f2(t). The timescale is 
t ~ To(Mfl)~ 6 , where a dominant iw-mode instability 
is assumed [61|. The mass can be determined as follows. 
The energy carried by the gravitational wave is (69j 



d 2 E 
dtdtt 



lim 



D 2 
16tt 



(h\ + hl) 



(4.8) 



Assuming that all the energy carried by the gravita- 
tional wave is extracted from the star, it follows dE/dt = 
—dM/dt. The angular momentum radiated in the az- 
imuthal mode m can be obtained from dE = LodJ/m. 
If this angular momentum is also completely extracted 
from the star, then J = IQ, where the moment of inertia 
/ = 2f3M(t) 3 is approximately constant in time. Setting 
dJ/dt — m/cudE/dt = 2/ (aft)dM/dt, the mass variation 
rate is 



M 
M 3 



Integration of Eq. (I4.9[) yields 



M~ 2 — M~ 2 

n 2 = nl + M ° f 



(4.9) 



(4-10) 



The mass and the angular velocity can be obtained 
by solving the previous equations with the condition 
dE/dt — —dM/dt. A typical solution is shown in Figure 
1 



V. DISCUSSION 

The above results show that ultra-compact objects 
with high redshift at their surface are unstable when 
rapidly spinning. This strengthens the role of BHs as 
candidates for astrophysical observations of rapidly spin- 
ning compact objects. 

Boson stars and gravastars easily develop ergoregion 
instabilities. Analytical and numerical results indicate 
that these objects are unstable against scalar field per- 
turbations. Their instability timescale is many orders of 
magnitude stronger than the instability timescale for or- 
dinary stars with uniform density [52] . In the large I = m 
approximation, suitable for a WKB treatment, gravita- 
tional and scalar perturbations have similar instability 
timescales. In the low-m regime gravitational perturba- 
tions are expected to have shorter instability timescales 
than scalar perturbations. 

For J > QAM 2 instability timescales can be as low 
as a few tenths of a second for solar mass objects and 
about a week for supermassive BHs, monotonically de- 
creasing for larger rotations. Therefore, high rotation is 
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evolve through the ergoregion instability, gravitational- 
wave detectors could easily identify them with a matched- 
filtering search. 

It would be interesting to repeat the above analy- 
sis for other ultra-compact objects, such as wormholes 
[72L l73l l74l . [75l . superspinars [76| and quark or fermion- 
boson stars [771 ] . The general arguments presented above 
suggest that these objects should be unstable. It would 
also be interesting to perform a more detailed analysis 
of the gravitational sector of the ergoregion instability 
and derive more refined templates for matched-filtering 
searches. A first step in this direction would be to com- 
pute axial perturbations for any value of the azimuthal 
number m along the lines of previous works [6l[ . 



t/M 
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FIG. 8: Evolution of mass, angular velocity and waveform 
of a thin-shell gravastar with ho = 10~ 4 , MQo = 0.25 and 
t/M = 0.041 during the linear perturbation phase. The mass 
difference is rescaled by a factor 10 for better visualization. 

an indirect evidence for horizons. The spin of an astro- 
physical compact object can be estimated by looking at 
the gas accreting near its surface 0, [7l| . A handful of 
fast-spinning BH-like objects have been reported @, 0]. 
The results of this paper suggest that these objects must 
indeed be BHs. 

The ergoregion instability evolution is characterized by 
very distinct waveforms. If compact astrophysical objects 
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